Effect of Temperature Wave on the Gas Transport in Liquid-Saturated Porous Media 
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We study the effect of surface temperature oscillations on gas mass transport through liquid- 
saturated porous media. Temperature wave induced by these oscillations and decaying deep in the 
massif creates the gas solubility wave along with the corresponding solute diffusion flux wave. When 
bubbles are immobilized by the surface tension force the only remaining mechanisms of gas mass 
transport are related to solute flux through liquid in pores. We evaluate analytically the generated 
time-average mass flux for the case of medium everywhere littered with gas bubbles and reveal 
the significant effect of the temperature wave on the gas release from the massif and bubble mass 
redistribution within the massif. Analytical theory is validated with numerical calculations. 

PACS numbers: 47.55.db, 66.10.C-, 92.40.Kf 



I. INTRODUCTION 

Transport of gases, as well as of any weakly soluble 
substances, in liquid-saturated porous media appears to 
possess unique features tl|-i3i|- Specifically, for a pore di- 
ameter small enough the pore-sized bubbles are immo- 
bilized by the surface tension force. The critical pore 
size can be readily accessed from the balance of the sur- 
face tension and buoyancy forces, which yields the value 
of order of 1 mm for the air- water system. Meanwhile, 
the bubbles which are large compared to the pore size 
will experience splitting during displacement due to the 
fingering instability 0|- As a result, the hydrodynamic 
transport of the gas phase through porous media becomes 
practically impossible for a small volumetric fraction of 
gas in pores (or gas saturation). The critical immobili- 
sation gas saturation varies from system to system and 
depends on the time scale (e.g., the leakage process sig- 
nificant on the time scales of millions of years is negligible 
for processes evolving in days), but remains within the 
range from 0.5-1% Q to several percent (see, e.g., 
The critical immobilized gas mass is at least by one or- 
der of magnitude larger than the feasible variations of 
gas mass dissolved in a saturated aqueous solution of 
oxygen, nitrogen, or methane under the Earth's near- 
surface conditions for pressure and temperature. The 
latter conditions are relevant for the processes in wet- 
lands and peatbogs. Thus the principal mechanism of 
the gas mass transport in groundwater turns out to be 
the transport of gas molecules dissolved in water. 

In the massif zone where the aqueous solution is 
saturated — i.e., porous massif is littered with bubbles 
of gas phase — the gas flux is solely determined by the 
field of solubility, which depends on pressure and tem- 
perature [1, 01 • In this relation it is remarkable that the 
temperature increase from 0°C to 20° C leads to the de- 
crease of the solubility of the main atmospheric gases and 
methane by factor 1.5. This means that the annual tem- 
perature wave propagating into the porous massif (e.g., 
see [1] ) gives rise to the solubility wave of significant am- 
plitude and associated diffusive fluxes. Hence, the annual 



temperature wave can signiflcantly affect the processes 
of saturating groundwater with gases or methane release 
from wetlands and peatbogs. 

The phenomenon under consideration is common and 
can be significant for various systems with different ori- 
gins of the surface temperature oscillations, including 
technological systems (filters, porous bodies of nuclear 
and chemical reactors, etc.). However, for the sake of 
convenience, in this paper, we first focus on the case fea- 
tured by the hydrostatic pressure gradient which is sig- 
nificant for geological systems, where pressure doubles 
on the depth of 10 meter, leading to significant change of 
solubility. The no pressure gradient case can be derived 
from the case of hydrostatic pressure. 

We restrict our consideration to the case when the so- 
lution is saturated everywhere in the massif. This re- 
striction is suggested not only by the motive of funda- 
mental theoretical interest related to problem novelty 
but also by its practical relevance. It is relevant for the 
methane release from peatbogs, as the process of methane 
generation there can be intense enough to maintain the 
saturation Q, and for atmospheric gases. For the lat- 
ter, it is important that the present time is the warm 
epoch against the background of the "main" cold state 
of 100 000-year Glacial-Interglacial cycles [l3|- During 
the cold period the groundwater is saturated with atmo- 
spheric gases under enhanced solubility conditions (low 
temperature), while during the warm period of dimin- 
ished solubility the excessive gas mass from the solution 
forms gaseous bubbles in pores. For low-permeability 
fractureless massifs, the dominant transport mechanism 
is the molecular diffusion. Molecular diffusion operates 
with equal efficiency downwards, during the cold periods 
when the groundwater is being saturated with gases, and 
upwards, during the warm periods when the system is 
losing the gas mass. The longer the saturating period is, 
the longer duration of the release period is needed to re- 
move the excessive amount of gases. Due to asymmetry 
between cold (long) and warm (short) periods groundwa- 
ter in the near-surface layers should remain gas-saturated 
during the warm period. Thus, for the current time, one 
should observe gradual release of the previously saturated 
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FIG. 1: (Color online) Average temperatures in New York 
City (circles) and their harmonic approximation (solid line) 



gases from massifs under the saturation condition. Short- 
time temperature waves can affect the rate of this release. 
Thus, we can see the relevance of the problem statement 
assuming the solution to be everywhere saturated. 

In this paper, we calculate the diffusion flux of gas in 
the presence of temperature wave and redistribution of 
the gas in non-dissolved gaseous phase. Both the cases 
of molecular diffusion and hydrodynamic dispersion are 
considered. The analytical results are validated with nu- 
merical simulation. The periods of negative tempera- 
tures with frozen groundwater are beyond the scope of 
this study and will be considered elsewhere. 



II. TEMPERATURE WAVE IN BUBBLY MEDIA 

We adopt harmonic annual oscillation of the surface 
temperature, Tq + Qq cos ujt, where Tq is the annual- 
mean temperature, Qq is the oscillation amplitude, and 
uj — 27r/year. This assumption is accurate enough due 
to two reasons: (i) real temperature records slightly de- 
viate from their harmonic reduction (see Fig.[T]) and (ii) 
the penetration depth of higher harmonics into massif 
rapidly decreases (see Eq. ^ below). With the specified 
surface temperature the temperature field in the half- 
space, which is governed by the heat diffusion equation 
dT/dt = with no heat flux condition at infinity, is 



T(z) = To + 606"'^'= cos(wt - /cz) , k=yU7/2x. (f) 

Here x is the heat diffusivity, the z-axis is oriented down- 
wards and its origin is on the massif surface. The pressure 
field is hydrostatic 
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where Pq is the atmospheric pressure, p is the 
(mostly, water) density, and g is the gravity. 

For pressure up to f atm and far from the solvent boil- 
ing temperature the solubility depends on temperature T 
and pressure P as follows 
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TABLE I: Chemical physical properties of nitrogen, oxygen, 
methane, and carbon dioxide (see text for details) . Eq. Q 
with q and X'"' (To, Po) specified in the table matches the 
experimental data from [1^ . Eq. ([5]) with Rd and v specified 
matches the experimental data from [ij. 







O2 


CH4 


CO2 


q = -G./fcB (K) 


781 


831 


1138 


1850 


X(°)(20°C, latin) (10"^) 


1.20 


2.41 


2.60 


68.7 


Ri (10-^° m) 


1.48 


1.29 


1.91 


1.57 


u (10"^Pa-s) 


9.79 


16.3 


28.3 


4.68 



where molar solubility X^'^^ is the molar amount of so- 
lute per 1 mole of solvent, Tq and Pq are reference values, 
the choice of which is guided merely by convenience rea- 
son, and X^'^\Tq,Pq) is the solubility at the reference 
temperature and pressure; the parameter q = —Gi/k-Q, 
with Gi being the interaction energy between a solute 
molecule and the surrounding solvent molecules, is pro- 
vided in Table m 

The instantaneous flux of molar fraction X in bulk 
obeys 
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where a is the thermodiffusion constant [llj . 

Generally, not only solubility but all material proper- 
ties of the system depend on temperature and pressure, 
flowever, feasible amplitudes of relative variations of the 
absolute temperature are below 5%. Hence, one can ne- 
glect variation of those characteristics which depend on 
temperature polynomially and consider variation of only 
those characteristics which depend on temperature expo- 
nentially: solubility ([3]) and the molecular diffusion co- 
efficient D. Further, the only characteristic sensitive to 
pressure below hundreds of atmospheres is gas solubility. 

To complete the mathematical formulation of the prob- 
lem, one needs to specify the dependence of molecular 
diffusion on temperature [IH 



D = 



(5) 



where fee is the Boltzmann constant, fj, is the dynamic 
viscosity of the solvent, Rd is the effective radius of the 
solute molecules with the "coefficient of sliding friction" 
/3, 1/ = RdP/i. The dependence of dynamic viscosity 
on temperature can be described by a modified Frenkel 
formula \\% 



M = Mo exp 



T 



(6) 



(3) 



For water, coefficient /Ltp = 2.42 -10 ^ Pa • s, a = W/k^ 
570 K {W is activation energy), and r = — 140K. 



3 



III. AVERAGE GAS MASS TRANSPORT 

Since the diffusion transport in liquids is by several or- 
ders of magnitude slower than the heat transfer, it can 
be well described in terms of average values over the tem- 
perature oscillation period. With the solubility law (|3|) 
diffusion flux ^ reads 



With expansion 
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D = Do + Di{T- To) + -D2{T - Tof 



(7) 



where Dq, Di, and D2 can be plainly evaluated from 
Eqs. ([5]) and straightforward but laborious analytical 
calculations yield mean flux 
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Beyond the penetration zone of the decaying temper- 
ature wave, i.e., for kz > 2, the diffusion flux ^ is ho- 
mogeneous, no gas mass accumulation or substraction 
is created by the flux divergence. Meantime, in the 
near-surface zone where the temperature wave is non- 
small, the flux divergence is non-zero and one observes 
the growth (or dissolution) of the gaseous phase; 



d{X,) d{J) 



dt 



dz 



(9) 



Here Xf, is the molar amount of matter in the gaseous 
phase per 1 mole of all the matter in pores. The latter 
equation is valid for small volumetric fraction of bubbles 
in pores, which well corresponds to the systems we con- 
sider. 

In Fig. [2] one can see that analytical results involving 
the expansion ([7]) match well the results of numerical 
calculations even for surface temperature oscillation am- 
plitude Qq as large as 15 K. Henceforth, our treatment 
relies entirely on these analytical results. 



IV. MOLECULAR DIFFUSION AND 
HYDRODYNAMIC DISPERSION 

Up to this point our results are derived for the molecu- 
lar diffusion mechanism. With molecular diffusion, terms 




To= 288 K, ©0= 15 K 



numeric 



— — — analytic 



12 3 4 

Depth below surface, z [m] 

FIG. 2: (Color online) Comparison of the numerically eval- 
uated divergence of mean flux with diffusion coefficient ([5]) 
and analytical expressions (|8]) and ([9| for annual tempera- 
ture wave 



{DiTqJDq) and {D2TQ / Dq) are of the order of magnitude 
of 10 [iTj; they make principle contribution to constants 
Ai and A2, and play a decisive role in the systems evo- 
lution. Noticeably, calculations reveal the effect of ther- 
modiffusion (a) to be practically unobservable against 
the background of other contributions to the diffusion 
flux. Hence, the uncertainty of the value of a for aque- 
ous gas solutions [l^ is not a significant issue for the 
results provided in this paper. 

However, along with the molecular-diffusion domi- 
nated systems, there are systems where horizontal filtra- 
tion flux of groundwater is present. This flux is treated 
as horizontal because it is parallel to the surface and 
natural systems are typically much more uniform in the 
horizontal directions than in the vertical one. Due to 
the microscopic irregularity of pore geometry the filtra- 
tion flux does fluid mixing which operates like an addi- 
tional diffusion mechanism, the "hydrodynamic disper- 
sion" [l|, [1^, [1^ . Although hydrodynamic dispersion was 
previously studied in the relation to the evolution of at- 
mosphere gases bubbles in aquifers [H , the effect of tem- 
perature waves was not addressed. 

When hydrodynamic dispersion is present, it is typi- 
cally by several orders of magnitude stronger than molec- 
ular diffusion (P, [l^, cf. in Fig. [3]), and the latter can 
be neglected. There is no analog of thermodiffusion for 
hydrodynamic dispersion and a = for this case. Fur- 
thermore, hydrodynamic dispersion depends on the flux 
strength and pore geometry but not temperature, and 
the correlations between instantaneous flux oscillations 
and temperature are not obvious; therefore, D — Dq and 
Di = D2 — 0. Finally, for hydrodynamic-dispersion dom- 
inated system, Eqs. ([8]) and ([9]) yield 
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In this paper we analyze both molecular-diffusion dom- 
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FIG. 3: (Color online) The rate of bubble growth from saturated aqueous solution of methane subject to annual temperature 
wave. The transport mechanism is the molecular diffusion under annual-mean temperature To — 288 K (a) and To — 298 K 
(b), and the hydrodynamic dispersion under annual-mean temperature To = 288 K (c) and Tq = 298 K (d). The surface 
temperature oscillation amplitude Oq is specified in plots. For molecular diffusion, rough theoretical assessment a — 1.8 [3] 
is adopted, although parametric study reveals the result to be nearly insensitive to the specific realistic value of a. For 
hydrodynamic dispersion Dvort ~ 2- 10~^ m^/s, which typically corresponds, e.g., to the horizontal filtration flux of groundwater 
f aq = 0.02 cm/s through the medium with the pore size ~ 1 mm (cf | 



inated systems and hydrodynamic-dispersion dominated 
systems. 



V. RESULTS AND DISCUSSION 

In Fig. [3] one can see the results of calculations of 
the methane bubble growth rate in molecular-diffusion 
and hydrodynamic-dispersion dominated systems. These 
plots demonstrate several unique features of these cases. 
Let us discuss these features and their origins in detail. 

Noteworthy, with molecular diffusion of methane, one 
can observe both bubble dissolution for lower mean tem- 
peratures (Fig.[3K) and bubble growth for higher mean 
temperatures (Fig.[3]D). This is dominantly controlled by 
parameters {DiTq/ Dq) and {D2TQ / Dq) and may vary 
from gas to gas. Analytical expression (|8]) contains both 
exponential and linear in z terms and multipliers. With 
purely exponential terms, the shape of the bubble growth 
rate profile would be determined only by k [see Eq. dH)]. 
Owing to the linear in z multiplier in Eq. (|5]), this toler- 
ance of the profile shape to parameters Ai and A2 is bro- 
ken by the competition between terms e~^'^^ and ze~'^^^ . 



In particular. Figs. [3^ and[5)D demonstrate that the pene- 
tration depth of the bubble depletion zone is nearly twice 
bigger than the penetration depth of the growth zone. 

On the contrary, for hydrodynamic-dispersion domi- 
nated systems, only growth of bubbles is possible. Ac- 
cording to Eq. ([TUl) , the shape of the bubble growth rate 
profile is determined by the second z-derivative term. 
This derivative profile depends on fc = yju:/2x- The sec- 
ond z-derivative is everywhere positive for k > pg/Pg = 
0.1 m~^, i.e., for temperature wave penetration depth be- 
low 10 m, while for k < pg/Po it turns negative in the 
upper zone above = — Po/pg. Since for natural 
systems the penetration depth of the annual temperature 
wave never reaches 10 m |8l|, the latter case, with bubble 
depletion above and growth below z*, is not feasible 
and can be relevant only for longer climate cycles. 

We discussed our analytical results in relation to geo- 
logical systems featured by the hydrostatic pressure gra- 
dient and annual temperature oscillations. Meanwhile, 
Eqs. ([8]) and (|10p hold valid for the systems where the 
hydrostatic pressure gradient is either absent or insignif- 
icant. The latter case may be, for example, daily temper- 
ature oscillation, the penetration depth of which is 0.1m 
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for soils — a depth on which the hydrostatic pressure in- 
crease is negligible. The case of no pressure gradient can 
be derived from Eq. ([8]) by eliminating the gravity; one 
finds 

(J)^iJoX(")(To,Po)A /""^^, +o(^^j. (11) 

Here we can see a non-trivial temperature wave effect on 
gas release for various technological and natural systems 
as well. 



VI. CONCLUSION 

We have addressed the problem of influence of the sur- 
face temperature oscillation (and consequent tempera- 
ture wave) on the gas transport through liquid-saturated 
porous media. Specifically, we have considered the prob- 
lem for the case of saturated gas solution, i.e., when 
bubbles are spread everywhere in the medium. Strong 
exponential dependence of the solubility and the molec- 
ular diffusion coefficient on temperature has been found 
to lead to non-negligible instantaneous diffusion fluxes 
although the amplitude of the relative variation of abso- 
lute temperature does not exceed 5%. Due to nonlinear- 



ity these fluxes are not averaged-out to zero but create 
the mean mass flux. 

Our main findings are expressed by Eq. and 
its reductions for hydrodynamic-dispersion dominated 
systems [Eq. ([TO)) ] and no pressure gradient systems 
[Eq. (ITT]) ]. These analytical equations have been shown 
to fairly match the results of numerical calculations (see 
Fig. 12). We have revealed that the temperature wave can 
either enhance or deplete the near-surface bubbly zone 
in molecular-diffusion dominated systems and only en- 
hance this zone for hydrodynamic-dispersion dominated 
systems (for instance, see Fig.jS]). 

The phenomenon we have addressed is expected to be 
of significance not only for gases but also for any weakly 
soluble substance, given its solubility is sensitive to tem- 
perature. 
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